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The emission of gravitational waves from a system of massive objects interacting on elliptical, 
hyperbolic and parabolic orbits is studied in the quadrupole approximation. Analytical expressions 
are then derived for the gravitational wave luminosity, the total energy output and gravitational 
radiation amplitude. A crude estimate of the expected number of events towards peculiar targets 
(i.e. globular clusters) is also given. In particular, the rate of events per year is obtained for the 
QO . dense stellar cluster at the Galactic Center. 
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I. INTRODUCTION 

m 

Gravitational-wave (GW) science has entered a new era. Experimentally, several GW ground-based-laser- 
interferometer detectors (10~ 1 kHz) have been built in the United States (LIGO) Europe (VIRGO and GEO) 
and Japan (TAMA) g|, and are now taking data at design sensitivity. Advanced optical configurations capable 
O r of reaching sensitivities slightly above and even below the so-called standard-quantum-limit for a free test-particle, 
have been designed for second [fj|] and third generation [6| GW detectors. A laser-interferometer space antenna (LISA) 
(1(T 4 - W- 2 Hz) mi ght fly within the next decade. Resonant-bar detectors (~ lfci?z)[8( are improving more and 
more their sensitivity, broadening their frequency band. At much lower frequencies, ~ 10 -17 Hz, future cosmic mi- 
crowave background (CMB) probes are devoted to detect GWs by measuring the CMB polarization Millisecond 
pulsar timing can set interesting upper limits [Ic| in the frequency range 10 -9 ~ 10~ s Hz. In this frequency range, 
j the large number of millisecond pulsars which will be detectable with the square kilometer array [ll| , would provide 
• an ensemble of clocks that can be used as multiple arms of a GW detector. 

' From a theoretical point of view, recent years have been characterized by numerous major advances due, essentially, 
to the development of numerical gravity. Concerning the most promising sources to be detected, the GW generation 
problem has improved significantly in relation to the dynamics of binary and multiple systems of compact objects as 
neutron stars and black holes. 

^,0 • Besides, the problem of non-geodesic motion of particles in curved spacetime has been developed considering the 
emission of GWs [l2l[l3l]. Solving this problem is of considerable importance in order to predict the accurate waveforms 
of GWs emitted by extreme mass-ratio binaries, which are among the most promising sources for LISA [l4j]. To this 
aim, searching for criteria to classify the ways in which sources collide is of fundamental importance. A first rough 
, criterion can be the classification of stellar encounters in collisional as in the globular clusters and in collisionless as 
in the galaxies [25| . A fundamental parameter is the richness and the density of the stellar system and so, obviously, 
we expect a large production of GWs in rich and dense systems. 

Systems with these features are the globular clusters and the galaxy centers. In particular, one can take into account 
the stars (early-type and late-type) which are around our Galactic Center, e.g. Sagittarius A* (SgrA*) which could 
be very interesting targets for the above mentioned ground-based and space-based detectors. 

In recent years, detailed information has been achieved for kinematics and dynamics of stars moving in the gravita- 
tional field of such a central object. The statistical properties of spatial and kincmatical distributions are of particular 
interest (see e.g. [H, [3 H3| ) ■ Using them, it is possible to give a quite accurate estimate of the mass and the size of the 
central object: we have (2.61 ± 0.76) x 10 6 M Q ) concentrated within a radius of 0.016pc (about 30 light-davs) flU [l9| . 
More precisely, in [18| , it is described a campaign of observations where velocity measurements in the central arcsec 2 
are extremely accurate. Then from this bulk of data, considering a field of resolved stars whose proper motions are 
accurately known, one can classify orbital motions and deduce, in principle, the rate of production of GWs according 
to the different types of orbits. This motivates this paper in which, by a classification of orbits in accordance with the 
conditions of motion, we want to calculate the GW luminosity for the different types of stellar encounters. A similar 
approach has been developed in [201 ] but, in that case, only hyperbolic trajectories have been considered. 

Following the method outlined in [2l|, [22j , we investigate the GW emission by binary systems in the quadrupole 
approximation considering bounded (circular or elliptical) and unbounded (parabolic or hyperbolic) orbits. Obviously, 
the main parameter is the approaching energy of the stars in the system (see also [23[ and references therein). We 
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expect that gravitational waves are emitted with a "peculiar" signature related to the encounter- type: such a signature 
has to be a "burst" wave- form with a maximum in correspondence of the periastron distance. The problem of 
bremsstrahlung-\ike gravitational wave emission has been studied in detail by Kovacs and Thorne [24] by considering 
stars interacting on unbounded and bounded orbits. In this paper, we face this problem discussing in detail the 
dynamics of such a phenomenon which could greatly improve the statistics of possible GW sources. 

The paper is organized as follows: in Sec. II, the main features of stellar encounters and orbit classification are 
reviewed. Sec. Ill is devoted to the emission and luminosity of GWs from binary systems in the different kinds of 
orbits assuming the quadrupole approximation. A discussion of the wave-form dependence from the orbital parameters 
is given in Sec. IV. In Sec. V, we derive the expected rate of events assuming the Galactic Center as a target. Section 
VI is devoted to concluding remarks. 



II. ORBITS IN STELLAR ENCOUNTERS 



Let us take into account the Newtonian theory of orbits since stellar systems, also if at high densities and constituted 
by compact objects, can be usually assumed in Newtonian regime. We give here a self-contained summary of the 
well-known orbital types in order to achieve below a clear classification of the possible GW emissions. We refer to the 
text books [25|, [2?| for a detailed discussion. 

A mass mi is moving in the gravitational potential $ generated by a second mass m,i. The vector radius and the 
polar angle depend on time as a consequence of the star motion, i.e. r = r(t) and <fi — 4>(t). With this choice, the 
velocity v of the mass mi can be parameterized as 

v = Vrr + v^cj) , (1) 
where the radial and tangent components of the velocity are, respectively, 

dr dej) 

In this case, the total energy and the angular momentum, read out 

1 (dr\ 2 L 2 7 

and 

^% M 

respectively, where fi = ^p^f- is the reduced mass of the system and 7 = Gmim 2 - 

We can split the kinetic energy into two terms where, due to the conservation of angular momentum, the second 
one is a function of r only. An effective potential energy V e //, 

Veff = ^ - 1 (5) 

is immediately defined. The first term corresponds to a repulsive force, called the angular momentum barrier. The 
second term is the gravitational attraction. The interplay between attraction and repulsion is such that the effective 
potential energy has a minimum. Indeed, differentiating with respect to r one finds that the minimum lies at ro = — 
and that 

2 
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Therefore, since the radial part of kinetic energy, 
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Kr = - n ,[ Jt ) (7) 



is non-negative, the total energy must be not less than VJj™, i.e. 
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The equal sign corresponds to the radial motion. For E m i n < E < 0, the trajectory lies between a smallest value r n 
and greatest value r max which can be found from the condition E = V e ft, i.e. 



7 / ( 1 \ 2 -^ 2 
r {mtn ,ma X } = ~2^± y [l^j) +^ ( 9 ) 

where the upper (lower) sign corresponds to r max (r min ). Only for E > 0, the upper sign gives an acceptable value; 
the second root is negative and must be rejected. 

Let us now proceed in solving the differential equations ((3|) and (j4]). We have 



dr dr dcj> L dr L d f 1 
dt d(f) dt [ir 2 d(f> (i d(f) \r 

and defining, as standard, the auxiliary variable u — 1/r, Eq. ([3]) takes the form 



(10) 



where u' = du/d<f) and we have divided by L 2 /2fi. Differentiating with respect to 4>, we get 

u'(u" + u-2£)=0 (12) 
hence either v! — 0, corresponding to the circular motion, or 



(13) 



which has the solution 



or, reverting the variable, 



u= ^r + Ccos(cj) + a) (14) 



2^ + Ccos((j) + a)\ (15) 

which is the canonical form of conic sections in polar coordinates [26[ ■ The constant C and a are two integration 
constants of the second order differential equation (fT3"|) . The solution lfT5"|) must satisfy the first order differential 
equation (fTTj) . Substituting (fTS"]) into (jTTJ) we find, after a little algebra, 

and therefore, taking account of Eq. ([5]), we get C 2 > 0. This implies the four kinds of orbits given in Table I 



° = ir + m (16) 



A. Circular Orbits 



Circular motion corresponds to the condition u' = by which one find r = L 2 / ^ where V e ff has its minimum. 
We also note that the expression for ro together with Eq.® gives 

ro = -^~ (17) 



Thus the two bodies move in concentric circles with radii, inversely proportional to their masses and are always in 
opposition. 
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c = o 


E — E m i n 


circular orbits 


o< \c\ < ft 


Emin < E < 


elliptic orbits 


\c\ = ft 


E = 


parabolic orbits 


\C\ > ft 


£ > 0, 


hyperbolic orbits 



Table I: Orbits in Newtonian regime classified by the approaching energy. 



B. Elliptical Orbits 

For < \C\ < fij/L 2 , r remains finite for all values of <f>. Since r{4> + 2ir) = r((f>), the trajectory is closed and it is 
an ellipse. If one chooses a = 0, the major axis of the ellipse corresponds to </> = 0. We get 



H*=o 



7M 
L 2 



C 



and 



(18) 



ru- 



7r ' max 



n 



and since r 



max i ' mm 



/■,,,,,, = 2a, where a is the semi-major axis of the ellipse, one obtains 

-l 



7/i 



" ^*| Train jj^ 

C can be eliminated from the latter equation and Eq.( [T5|) and then 



(19) 



_L 

2E 



(20) 



Furthermore, if we denote the distance r^ = ^^ 2 by h the so-called semi-lotus rectum ox the parameter of the ellipse, 
we get 



7/i 



(21) 



and hence the equation of the trajectory 



1 + e cos < 



(22) 



where e = \ ^—^ is the eccentricity of the ellipse. 



C. Parabolic and Hyperbolic Orbits 



These solutions can be dealt together. They correspond to E > which is the condition to obtain unbounded 
orbits. Equivalently, one has |C| > 7/i/L 2 . 
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The trajectory is 

r = /(l + ecos0) _1 (23) 

where e > 1. The equal sign corresponds to E — . Therefore, in order to ensure positivity of r, the polar angle <p 
has to be restricted to the range given by 

l + ecos0>O (24) 

This means cos0 > — 1, i.e. <f> € (— and the trajectory is not closed any more. For (f> — > ±tt, we have r — ► oo. 
The curve (|23|) . with e = 1, is a parabola. For e > 1, the allowed interval of polar angles is smaller than <fi € (— 7r, 7r) , 
and the trajectory is a hyperbola. Such trajectories correspond to non-returning objects. 

III. GRAVITATIONAL WAVE LUMINOSITY IN THE QUADRUPOLE APPROXIMATION 

At this point, considering the orbit equations, we want to classify the gravitational radiation for the different stellar 
encounters. To this aim, let us start with a short review of the quadrupole approximation for GW radiation. We add 
this discussion for the sake of completeness, but we send the Reader to the References [28|, HiJ HO, HH for a detailed 
exposition. 

The Einstein field equations give a description of how the curvature of space-time is related to the energy-momentum 
distribution. In the weak field approximation, moving massive objects produce gravitational waves which propagate 
in the vacuum with the speed of light. In this approximation, we have 

9nv = Spiv + Kh^, (25) 
whit k the gravitational coupling. The field equations are 



□V, = ~o kT ^ ( 26 ) 



where 



h>n V — — -S^hxx , (27) 

and Tuv is the total stress-momentum-energy tensor of the source, including the gravitational stresses. 
A plane GW can be written as 

h^u = hfj, v = he^v cos(ujt - k • x) (28) 

where h is the amplitude, lj the frequency, k the wave number and e^ t „ is a unit polarization tensor, obeying the 
conditions 

Let us assume a gauge in which e^„ is space-like and transverse; thus, a wave travelling in the z direction has two 
possible independent polarizations: 

ei = -j=(xx - yy) e 2 = -j=(xy - yx). (30) 

One can now search for wave solutions of |26|) from a system of masses undergoing arbitrary motions, and then 
obtain the power radiated. The result, assuming the source dimensions very small with respect to the wavelengths 

(quadrupole approximation [2jJ), is that the power —— radiated in a solid angle f2 with polarization is 
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dE 



( dPQi 



G 



8nc 5 V dt 3 



where Qij is the quadrupole mass tensor 



Qij — ^""^ m a(3x a x J a Sijr a ) , 



(31) 



(32) 



G being the Newton constant, r a the modulus of the vector radius of the a — th particle and the sum running over all 
masses m a in the system. It has to be noted that the result is independent of the kind of stresses which are present 
into the dynamics. If one sums (|3ip over the two allowed polarizations, one obtains 



E 

pol 



dE 

dn 



G 
8ttc 5 

1 



cPQij d 3 Q i: 
dt 3 dt 3 



d 3 Qij d 3 Q k j 
2n, . „ J n k 



dt 3 



d 3 Q t 
dt 3 



-TljTlk 



dt 3 

d 3 Q jk 
dt 3 



1 ( d 3 Q v 

2 I dt 3 



(33) 



where n is the unit vector in the radiation direction. The total radiation rate is obtained by integrating (|33[) over 
all directions of emission; the result is 



«_e «(Qg'Q' 3 " J ) .... 

dt 45c 5 { ' 

where the index (3) represents the differentiation with respect to time, the symbol <> indicates that the quantity 
is averaged over several wavelengths. This crucial point is linked with the difficulties of localizing gravitational energy 
so the right hand side of Eq. ([34]) cannot be viewed as an instantaneous quantity. This problem has been already 
faced for circular and elliptical orbits in [3, H2, [33|. For hyperbolic and parabolic orbits, it is crucial to estimate 
the quantity in the right hand side of Eq. (l34D in the zone where stars are slightly changing their trajectories, that 
means at peri-astron, while we expect no emission in asymptotic regime of stars approaching to or going away from 
this region. For a detailed discussion in the hyperbolic case, see (20j . 

With this formalism, it is possible to estimate the amount of energy emitted in the form of GWs from a system of 
massive objects interacting among them [2~H l22j . In this case, the components of the quadrupole mass tensor in the 
equatorial plane (9 = tt/2) are 

/xr 2 (3cos 2 — 1) , 
[ir 2 (3sm 2 <f>- 1) , 
-lir 2 , 

(35) 

Q zx — , 



Qxx 
Qyy 
Q zz 
Qxz 



Qyz — Qzy i 

Qxy = Qyx = cos (j) sin (/) , 

where the masses m\ and have polar coordinates \vi cos cos 0, r% cos#sin</), rising} with i = 1,2. The origin 
of the motions is taken at the center of mass. Such components can be differentiated with respect to the time as in 
Eq. (|34p . In doing so, we can use some useful relations derived in the previous Section. 
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A. GW luminosity from circular and elliptical orbits 

Using Eq: ([22)l . let us derive the angular velocity equation 



• JGl(mi + m 2 )(ecos(j) + l) 2 

0= p 

and then, from (|35[) . the quadrupolar components for the elliptical orbits 



(36) 



d 3 Q 

XX 

dt 3 

d Qyy 

dt 3 



d Qxy 

~di 3r 



= /3(24 cos + e(9 cos 20) + 1 1)) sin (37) 

-/3(24 cos + e(13 + 9 cos 2(f))) sin 0) (38) 

^£i = -2/? e sin0 (39) 

= /3(24 cos + e(ll + 9 cos 20)) sin 0) (40) 



where 



Being 



Gl(mi + m 2 )) 3 / 2 ^(e cos 0+1) 



2 



(41) 



Qg)Q(3)y = <J [ (mi + m2 )3 M 2 (1 + ecQS( 



Z 5 

(415e 2 + 3(8 cos + 3e cos 20) 
(72 cos + e(70 + 27 cos 20))) sin 2 ] 

the total power radiated is given by 

dE ° 3 

-dt = 45J^ m 

where 

/(0) = [(mi + m 2 )'V(l + ecos0) 4 (415e 2 + 3(8cos0 + 3ecos20) (72 cos + e(70 + 27cos 20))) sin 2 ] (42) 
The total energy emitted in the form of gravitational radiation, during the interaction, is given by 



AE 



POO 


dE 


Iq 


dt 



dt . (43) 



From Eq.Q, we can adopt the angle as a suitable integration variable. In this case, the energy emitted for 

01 < < 02 iS 



Afl " M ' j)= 4^' fU>)dn. (-1-1) 



and the total energy can be determined from the previous relation in the limits 0i — > and 02 — * n. Thus, one has 

AE = F{e) (45) 

where F{e) depends on the initial conditions only and is given by 

, N (13824 + 102448e 2 + 59412e 4 + 2549e 6 ) , 

F ^ = ~ 2880 1 (46) 

In other words, the gravitational wave luminosity strictly depends on the configuration and kinematics of the binary 
system. 
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B. GW luminosity from parabolic and hyperbolic orbits 



Also in this case, we use Eq: ([23|) and Eq. (|35j) to calculate the quadrupolar formula for parabolic and hyperbolic 
orbits. The angular velocity is 



and the derivative are 



Z 2 L(ecos0 + l) 2 (47) 



d 3 

= p(24cos0 + e(9cos20+ll))sin0 (48) 

dt i 



d Q yy 



= -p(24cos0 + e(13 + 9cos20))sin0) (49) 



dt 3 

d 3 Q 



dt 3 



= -2pesin0 (50) 



dt 3 2 



p(ecos0+l) (5ecos + 8 cos 20 + 3ecos30) (51) 



where 



p = I 4 L 3 p(e cos + l) 2 . (52) 

The radiated power is given by 

dE Gp 2 ([314e 2 + (1152cos(0 + 187ecos20 - 3(8Ocos30 + 3Oecos40 + 48 cos 50 + 9ecos60))e - 192 cos40 + 576] 
~dt ~ 120c 5 
then 

dE fflW 

~t = -^20^ m (53) 

where /(0), in this case, is 

/(0) = (314e 2 + (1 152 cos(0 + 187e cos 20 - 3(80 cos 30 + 30e cos 40 + 48 cos 50 + 9e cos 60))e - 192 cos 40 + 576) 

Then using Eq. (|34|) . the total energy emitted in the form of gravitational radiation during the interaction as a 
function of is given by 

, . G/ 8 i 6 7r (l271e 6 + 24276e 4 + 34768e 2 + 4608) p 2 , 
AE^t, 2 ) = ^ ^-g UL d( f> , (55) 

and the total energy can be determined from the previous relation in the limits 0i — > — tt and 02 — * 7r in the parabolic 
case. Thus, one has 

480c 5 w ' v 7 

where F(e) depends on the initial conditions only and is given by 

F(e) = (l271e 6 + 24276e 4 + 34768e 2 + 4608) . (57) 
In the hyperbolic case, we have that the total energy is determined in the limits 0! — ► — - — and 02 — > — — , i.e. 

r>l8 j-6.,2 
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where F(e) depends on the initial conditions only and is given by 

F(e) = [315tt (I271e 6 + 24276e 4 + 34768e 2 + 4608) + (59) 

+ 16e e e^926704%/2-7e(3319e 2 - 32632 V26 + 55200)) - 383460 + 352128V2 . (60) 

As above, the gravitational wave luminosity strictly depends on the configuration and kinematics of the binary 
system. 

IV. GRAVITATIONAL WAVE AMPLITUDE 

Direct signatures of gravitational radiation are its amplitude and its wave-form. In other words, the identification of 
a GW signal is strictly related to the accurate selection of the shape of wave-forms by interferometers or any possible 
detection tool. Such an achievement could give information on the nature of the GW source, on the propagating 
medium, and , in principle, on the gravitational theory producing such a radiation (35j . 

It is well known that the amplitude of GWs can be evaluated by 

1C 

h jk (t,R) = ^Q^ k , (61) 

R being the distance between the source and the observer and {j, k} — 1, 2. Let us now derive the GW amplitude in 
relation to the orbital shape of the binary systems. 

A. GW amplitude from elliptical orbits 

Considering a binary system and the single components of eq. (|61[) , it is straightforward to show that 

/ 11 - 2G G(mi+m 2 )Ai(13ecos0+12cos20+£(4e+3cos30)) 

;,22 _ 2G G(mi+m 2 )At(17ecos0+12cos20+e(8e+3cos30)) / £<i\ 

n ~ Rdf 21 ' W 

1 12 _ u21 _ 2G G(roi+ro 2 )At(13e cos 0+12 cos 20+e(4e+3 cos 30)) 
" — n ~ Rc 4 - 21 ' 

so that the expected strain amplitude h ~ (h 2 x + h 2 2 + ^h^) 1 ^ 2 turns out to be 

h = ° ^t^ 2 ^ (3(13e cos + 12 cos 20 + e(4e + 3 cos 30)) 2 + (17e cos 0+12 cos 20 + e(8e + 3 cos 30)) 2 ) i , (63) 

which, as before, strictly depends on the initial conditions of the stellar encounter. A remark is in order at this 
point. A monochromatic gravitational wave has, at most, two independent degrees of freedom. In fact, in the TT 
gauge, we have h+ = hxi + ft-22 and h x = hvz + ^21 (see e.g. [32]). As an example, the amplitude of gravitational 
wave is sketched in Fig. [1] for a stellar encounter close to the Galactic Center. The adopted initial parameters are 
typical of a close impact and are assumed to be b = 1 AU and vo = 200 Kms -1 , respectively. Here, we have fixed 
Mi = M 2 = 1.4Mq. The impact parameter is defined as L = bv where L is the angular momentum and v the 
incoming velocity. We have chosen a typical velocity of a star in the galaxy and we are considering, essentially, 
compact objects with masses comparable to the Chandrasekhar limit (~ 1.4M©). This choice is motivated by the 
fact that ground-based experiments like VIRGO or LIGO expect to detect typical GW emissions from the dynamics 
of these objects or from binary systems composed by them (see e.g. (30j|). 

B. GW amplitude from parabolic and hyperbolic orbits 

In this case the single components of Eq. (|6"Tj) for a parabolic and hyperbolic orbits, are 

h 11 = - G ^ 4 2/J (13e cos + 12 cos 20 + £ (4 £ + 3 cos 30)) , 

h 22 = S^Ji (I7 e cos + 12 cos 20 + e(8e + 3 cos 30)) , (64) 



h 12 = h 21 = - 3G ' i ^ (4 cos + e(cos 20 + 3)) sin 
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Figure 1: The gravitational wave-forms from elliptical orbits shown as function of the polar angle <j>. We have fixed Mi = M2 = 
IAMq. M2 is considered at rest while M\ is moving with initial velocity vo = 200 Kms" 1 and an impact parameter 6=1 AU. 
The distance of the GW source is assumed to be R = 8 kpc and the eccentricity is e = 0.2, 0.5, 0.7. 

and then the expected strain amplitude is 

2l 4 L 4 u 2 i 

h= 77f-(lOe 4 + 9e 3 cos30 + 59e 2 cos2(/) + 59e 2 + (47e 2 + lO8)ecos0 + 36) 2 , (65) 

c 4 ft 

which, as before, strictly depends on the initial conditions of the stellar encounter. We note that the gravitational 
wave amplitude has the same analytical expression for both cases and differs only for the value of e which is e = 1 
if the motion is parabolic and the polar angle range is 4> € (— 7r, 7r), while it is e > 1 and 4> 6 (— 7r,7r) for hyperbolic 
orbits. In these cases, we have non-returning objects. 

The amplitude of the gravitational wave is sketched in Figs. [2] and [3] for stellar encounters close to the Galactic 
Center. As above, we consider a close impact and assume 6=1 AU cm, vo = 200 Kms -1 and Mi = M2 = 1.4M Q . 

V. RATE AND EVENT NUMBER ESTIMATIONS 

An important remark is due at this point. A galaxy is a self-gravitating collisionless system where stellar impacts 
are very rare [25]. From the GW emission point of view, close orbital encounters, collisions and tidal interactions 
should be dealt on average if we want to investigate the gravitational radiation in a dense stellar system as we are 
going to do now. 

Let us give now an estimate of the stellar encounter rate producing GWs in some interesting astrophysical conditions 
like a typical globular cluster or towards the Galactic Center after we have discussed above the features distinguishing 
the various types of stellar encounters. Up to now, we have approximated stars as point masses. However, in dense 
regions of stellar systems, a star can pass so close to another that they raise tidal forces which dissipate their relative 
orbital kinetic energy. In some cases, the loss of energy can be so large that stars form binary or multiple systems; 
in other cases, the stars collide and coalesce into a single star; finally stars can exchange gravitational interaction in 
non-returning encounters. 

To investigate and parameterize all these effects, we have to compute the collision time t co ii, where l/t co u is the 
collision rate, that is, the average number of physical collisions that a given star suffers per unit time. For the sake 
of simplicity, we restrict to stellar clusters in which all stars have the same mass m. 

Let us consider an encounter with initial relative velocity Vo and impact parameter b. The angular momentum per 
unit mass of the reduced particle is L = bv a . At the distance of closest approach, which we denote by r co u, the radial 
velocity must be zero, and hence the angular momentum is L = r co uv max , where v max is the relative speed at r co u. 
From the energy equation ([3]), we have 
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Figure 2: The gravitational wave-forms for a parabolic encounter as a function of the polar angle <f>. As above, Mi = M2 = 
IAMq and M2 is considered at rest. Mi is moving with initial velocity Vo = 200 Kms" 1 with an impact parameter b — 1 AU. 
The distance of the GW source is assumed at R = 8 kpc. The eccentricity is e = 1. 




Figure 3: The gravitational wave-forms for hyperbolic encounters as function of the polar angle <f>. As above, we have fixed 
Mi = M2 = IAMq. M2 is considered at rest while Mi is moving with initial velocity vo = 200 Kms - and an impact parameter 
6=1 AU. The distance of the source is assumed at R — 8 kpc. The eccentricity is assumed with the values e = 1.2, 1.5, 1.7 . 



b 2 = 



4Gmr. 



coll 



coll 



(66) 



If we set r co u equal to the sum of the radii of two stars, then a collision will occur if and only if the impact parameter 
is less than the value of b, as determined by Eq. (f6"6"]) . 

Let /(v a )d 3 v a be the number of stars per unit volume with velocities in the range v Q + d 3 v a . The number of 
encounters per unit time with impact parameter less than b which are suffered by a given star is just /(v a )d 3 v a times 
the volume of the annulus with radius b and length vq, that is, 
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J /(v Q )7r6 2 « d 3 v Q (67) 



where vo = |v — v a | and v is the velocity of the considered star. The quantity in Eq. (|6"7| is equal to l/t co u for a 
star with velocity v: to obtain the mean value of l/t co u, we average over v by multiplying (|6T[) by f(v)/v, where 
v = j /(v)<i 3 v is the number density of stars and the integration is over d 3 v. Thus " 



(r coll |v - v.| + d 3 vd 3 v a (68) 



1 _ v_ f -t„?^?\ii^ ( , . AGmr, 

tcoll 87T 2 CT 



We now replace the variable v Q by V = v — v a . The argument of the exponential is then — 
and if we replace the variable v by v cm = v — — V (the center of mass velocity), then we have 

1 v f f„,2 ,t/2 W „ 2 / „ AGmr, 



fv _i V ) 2 + I^ 



e 



Uoii SttV 6 ./ V V 



(iC+v-W PooU y + ^ " co ' ; dV . (69) 



The integral over v cm is given by 



Thus 

1 TT 1 ^ 



t co /; 2cr 3 

The integrals can be easily calculated and then we find 



V /j ! j3„ ,^3/2 3 / 7n \ 

cm/ d V cm = 7T ' (7 . (7U) 





e -V 2 /4a 2 ( r 2 M y3 + AGmVrcM ) dV (71) 



= V7iwr co „ + . (72) 

tcoll a 

The first term of this result can be derived from the kinetic theory. The rate of interaction is vY> (V) , where £ is 
the cross-section and (V) is the mean relative speed. Substituting £ = 7rr 2 on and (V) = 4a/ ^tt (which is appropriate 
for a Maxwellian distribution whit dispersion a) we recover the first term of (|T2|) . The second term represents the 
enhancement in the collision rate by gravitational focusing, that is, the deflection of trajectories by the gravitational 
attraction of the two stars. 

If is the stellar radius, we may set r co ii = 27'*. It is convenient to introduce the escape speed from stellar surface, 

2Gm , . ^ ^ 
, and to rewrite hq. (|72J) as 

r = -!- = 16V^r 2 (l + ^r) = 16vWr 2 (l + 6), (73) 

tcoll V 4(7 / 

where 



v 



2 



Gm 



is the Safronov number (25|. In evaluating the rate, we are considering only those encounters producing gravitational 
waves, for example, in the LISA range, i.e. between 10 -4 and 10~ 2 Hz (see e.g. (36J). Numerically, we have 

i-io ( v \ ( a \ ( 10pc\ 3 . _i 



5.5 x lCT iu — - —— — — yrs- L 6 « 1 (75) 

VIOkms- 1 / \ UA 2 J V R ' 
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— °- 10 (To^) 2 (lo^Fr)(^)(T) 3 --' e». <*> 

If >> 1, the energy dissipated exceeds the relative kinetic energy of the colliding stars, and the stars coalesce 
into a single star. This new star may, in turn, collide and merge with other stars, thereby becoming very massive. As 
its mass increases, the collision time is shorten and then there may be runaway coalescence leading to the formation 
of a few supermassive objects per clusters. If << 1, much of the mass in the colliding stars may be liberated and 
forming new stars or a single supermassive objects (see [37l[38|). 

Note that when we have the effects of quasi-collisions in an encounter of two stars in which the minimum separation 
is several stellar radii, violent tides will raise on the surface of each star. The energy that excites the tides comes from 
the relative kinetic energy of the stars. This effect is important for >> 1 since the loss of small amount of kinetic 
energy may leave the two stars with negative total energy, that is, as a bounded binary system. Successive peri-center 
passages will dissipates more energy by GW radiation, until the binary orbit is nearly circular with a negligible or 
null GW radiation emission. 

Let us apply these considerations to the Galactic Center which can be modelled as a system of several compact 
stellar clusters, some of them similar to very compact globular clusters with high emission in X-rays (39j . 

For a typical compact stellar cluster around the Galactic Center, the expected event rate is of the order of 
2 x 10~ 9 yrs -1 which may be increased at least by a factor ~ 100 if one considers the number of globular clusters in 
the whole Galaxy eventually passing nearby the Galactic Center. If the compact stellar cluster at the Galactic Center 
is taken into account and assuming the total mass M ~ 3 x 10 6 M Q , the velocity dispersion a ~ 150 km s _1 and the 
radius of the object R ~ 10 pc (where = 4.3), one expects to have ~ 10 -5 open orbit encounters per year. On 
the other hand, if a cluster with total mass M ~ 10 6 M Q , a ~ 150 km s -1 and R ~ 0.1 pc is considered, an event 
rate number of the order of unity per year is obtained. These values could be realistically achieved by data coming 
from the forthcoming space interferometer LISA. As a secondary effect, the above wave-forms could constitute the 
"signature" to classify the different stellar encounters thanks to the differences of the shapes (see the above figures). 



VI. CONCLUDING REMARKS 



We have analyzed the gravitational wave emission coming from stellar encounters in Newtonian regime and in 
quadrupole approximation. In particular, we have taken into account the expected luminosity and the strain amplitude 
of gravitational radiation produced in tight impacts where two massive objects of IAMq closely interact at an impact 
distance of 1AU. Due to the high probability of such encounters inside rich stellar fields (e.g. globular clusters, bulges 
of galaxies and so on), the presented approach could highly contribute to enlarge the classes of gravitational wave 
sources (in particular, of dynamical phenomena capable of producing gravitational waves). In particular, a detailed 
theory of stellar orbits could improve the statistic of possible sources. 
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